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Abstract 

Utilizing the large dual description of a metastable system of branes and anti- 
branes wrapping rigid homologous S^'s in a non-compact Calabi-Yau threefold, we 
study phase transitions induced by changing the positions of the S'^'s. At leading 
order in 1/A^ the effective potential for this system is computed by the planar limit 
of an auxiliary matrix model. Beginning at the two loop correction, the degen- 
erate vacuum energy density of the discrete confining vacua split, and a potential 
is generated for the axion. Changing the relative positions of the S'^'s causes dis- 
crete jumps in the energetically preferred confining vacuum and can also obstruct 
direct brane/anti-brane annihilation processes. The branes must hop to nearby S'^'s 
before annihilating, thus significantly increasing the lifetime of the corresponding 
non-supersymmetric vacua. We also speculate that misaligned metastable glueball 
phases may generate a repulsive inter-brane force which stabilizes the radial mode 
present in compact Calabi-Yau threefolds. 



July, 2007 



*e-mail: jheckmaiiOf as .harvard, edu 
^e-mail: vaf aSphysics .harvard, edu 



1 Introduction 



Metastable vacua of supersymmetric string and field theories possess tfie attractive 
feature that in contrast to a generic non-supersymmetric system, the underlying 
supersymmetry of the theory often provides better control over the dynamics of 
the vacuum. String theory realizations of metastable vacua have been discussed 
in [1-11]. Recent progress in finding non-supersymmetric metastable vacua in su- 
persymmetric QCD-like field theories was achieved in [12] and subsequent string 
theory realizations of this work [13-27]. 

Combining these stringy insights with the dual closed string description of large 
open string systems such as the AdS / CFT correspondence [28-30] and geometric 
transitions [1,31,32], it was shown in [6] that much of the rigid structure oi M = 2 
supersymmetry remains intact for metastable vacua in type IIB string theory given by 
D5-branes and anti-D5-branes wrapping distinct minimal size S'^'s in a non-compact 
Calabi-Yau geometry of the form: 

y2 = W\xf + UV (1) 

where W'{x) = g{x — ai) ■ ■ ■ (x — a„) is a polynomial of degree n, the variables 
x,y,u,v E C and the minimal size S^'s are located at a; = a^. The Oj correspond to 
non-normalizable modes in the non-compact geometry which determine the relative 
separation between the branes. In the following we shall denote the number of 
branes (resp. anti-branes) by positive (resp. negative) integers. In the closed string 
dual description the branes (resp. anti-branes) wrapping S'^'s are replaced by S'^'s 
threaded by some amount of positive (resp. negative) RR flux. The geometry of 
the closed string dual description is given by blowing down the S"^ and introducing 
a complex deformation of equation ([1]) by a degree n — 1 polynomial in x: 

y^ = W'{x)^ + fn-l{x)+UV (2) 
n 

= g^Yl{x-ai){x-ar)+uv. (3) 

i=l 

This deformation splits the roots of W'{x)^ to and a~ for all i and thus defines a 
Riemann surface with n branch cuts located near each of these roots. 

At leading order in the 1/A^ expansion, the fluxes spontaneously break all of the 
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Figure 1: Depiction of flux line annihilation in a two cut geometry which initially 
consists of > units of flux through one cut and — through the other. In the 
depiction on the left, M > D5-branes wrap the interpolating 3-cycle between the 
two 5''^'s supported by flux. In the Minkowski spacetime the stack of D5-branes 
separates a bubble of vacuum with flux numbers (M — N,N — M) from one with 
flux numbers {—N,N). The end of the annihilation process and the corresponding 
flux numbers are shown on the right. 

M = 2 supersymmetry of the Calabi-Yau compactification. Due to the fact that the 
moduli space is still governed by rigid special geometry, the matrix model techniques 
developed in [33-36] still determine all higher order corrections to the form of the 
periods in the closed string dual. In the open string description these higher order 
corrections correspond to multi loop contributions to the glueball potential. In [8] 
a rich phase structure for the two cut system was uncovered by including two loop 
corrections to the effective potential. 

To frame some of the discussion to follow, we now review the two distinct ways 
in which the two cut system exhibits metastable behavior [8]. One decay mode is 
given by direct flux line annihilation where the domain wall delimiting the bubble 
of vacuum with lower flux is a stack of D5-branes wrapping the compact interpo- 
lating 3-cycle between the two S'^'s supported by flux. See figure [T] for a depiction 
of this process. Two loop effects reveal another way in which the vacua of the 
brane/anti-brane system are metastable. Whereas at one loop order the distinct 
confining vacua of the theory are energetically degenerate, beginning at two loop 
order this degeneracy is lifted. At large A^ and when the scale of confinement is not 
exponentially suppressed, the many confining vacua of the theory are all metastable. 
Geometrically, these vacua correspond to distinct orientations of the branch cuts. 
The decay to the lowest energy confining vacuum corresponds to a re-alignment of 
the branch cuts. See figure [2] for a depiction of this decay process. 

Although the analysis of brane/anti-brane configurations with two minimal size 
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Figure 2: Two loop corrections to the glueball potential lift the degeneracy in energy 
density between the confining vacua of the theory. These vacua correspond to 
distinct orientations of the branch cuts in the closed string dual geometry. The 
metastable branch cut orientations denoted by dashed lines eventually decay to the 
energetically preferred configuration denoted by a solid line. 

S^'s should provide an adequate description of "two body" interactions, a generic 
compact Calabi-Yau threefold will typically contain a large number of rigid homol- 
ogous S'^'s which may also be wrapped by branes and anti-branes. For example, 
the quintic Calabi-Yau threefold contains 2875 rational curves of precisely this type. 
Perhaps surprisingly, we find that the presence of additional unwrapped S'^'s intro- 
duces a further layer of phase structure. 

As already mentioned in the two cut case, just as a supersymmetric configuration 
of Ni branes contains Tr{—l)^ = Ni--- Nn energetically degenerate confining vacua, 
at one loop order the vacua of the brane/anti-brane system also remain energetically 
degenerate. It follows from an analysis similar to the two cut case that the two 
loop contribution to the energy density lifts this degeneracy. We also find, however, 
that changing the location of the minimal size S'^'s produces discrete jumps in the 
energetically preferred confining vacuum. 

Interpreting the two loop energy density as a potential for the axion of our theory, 
we next study the minima of this potential as a function of the Oj's. We find that the 
presence of additional unwrapped 5'^'s contributes a moduli space of values such that 
the effective ^-angle of each brane theory relaxes to zero at the minimum of the axion 
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potential. We next consider symmetric configurations where all S^'s are occupied 
and show that the discrete symmetries of the geometry translate into constraints on 
strong CP violation in brane/anti-brane configurations which remain invariant under 
the discrete symmetry x ^x. 

Although a general analysis of the critical points of the n-cut system appears 
quite difficult using present techniques, to keep the analysis tractable we study con- 
figurations with Ni = N D5-branes and N2 = —N anti-D5-branes wrapping minimal 
size S'^'s in a geometry defined by a collection of minimal size 5'^'s which all lie on the 
real axis and remain invariant under the map x ^ —x. We find that metastability 
is lost once the effective 't Hooft coupling becomes too large. At this point, the size 
of the corresponding branch cuts will begin to expand. A generic configuration of 
additional S^'s will typically obstruct a direct collision between the cuts supported 
by positive and negative flux. Such obstructions can lead to phase transitions to 
non-Kahler geometries of the type found in [8]. 

The additional S'^'s can also obstruct direct brane/anti-brane annihilation. Uti- 
lizing the partial classification of BPS states of geometrically engineered M = 2 
gauge theories obtained in [37], we show that the presence of additional minimal size 
S'^'s can sometimes cause the most efficient means of annihilation to proceed via a 
multi-domain wall process where a brane must first tunnel to an unoccupied minimal 
size S*^ before annihilating against an anti-brane. 

Combining the analysis of the previous sections, we also briefly comment on the 
stabilization of the radial mode which is present in compact Calabi-Yau threefolds. 
Whereas the one loop contribution to the vacuum energy density generates a loga- 
rithmic Coulomb attraction term, the two loop contribution generates a power law 
contribution which is repulsive for appropriate glueball phases. Although sugges- 
tive, we find that the value of the radial mode is very close to the value which would 
destabilize the size of the glueball fields. 

The rest of this paper is organized as follows. In section [2] we review the con- 
jectured geometric transition for the brane/anti-brane system and briefly discuss the 
phase structure found in [8]. In this same section we also include a short discussion 
on the regime of string theory parameters for which we expect to have an accurate 
description of the low energy dynamics. In section [3] we compute the two loop con- 
tribution to the period matrix of the n-cut geometry. We briefly discuss the form 
of the effective potential for geometries with unoccupied S'^'s in section HI We next 
compute corrections to the vacuum energy density as well as constraints on axion 



4 



potentials in sections [5] and [6], respectively. In section [7] we introduce a subclass of 
critical points in the presence of additional S^'s and study the breakdown of metasta- 
bility in section [HI In section [9] we study brane/anti-brane annihilation processes 
and in section [10] we speculate on the stabilization of the radial mode. Section [11] 
presents our conclusions. 

2 Geometrically Induced Metastability 

In this paper we consider metastable vacua of type IIB string theory. More pre- 
cisely, we study D5-branes and anti-D5-branes which fill Minkowski space and wrap 
n minimal size rigid homologous S'^'s of a local Calabi-Yau threefold defined by the 
hypersurface of equation ([T|). This configuration is metastable because the tension of 
the branes creates a potential barrier against moving off of the minimal size S"^. We 
now review some further properties of this system, closely following the discussion 
in [6,8]. At each point in the complex x-plane the area of the S*^ is: 

A{x)= (\W\x)\^ + \r\^y^ (4) 

where r denotes the stringy volume of the minimal size S'^'s given by turning on a 
non-trivial NS B2 field. The bare gauge coupling of the open string system is: 

47ri ^YM 47ri /• [ 

- a (Ao) = ^— + -— = B2+ C2 (5) 

9ym 27r g, J J 

where C2 denotes the RR two form. In the large limit, N 00 but Qs scales as 
A^~^ so that the associated 't Hooft coupling QsN remains finite. 

In the holographic dual description, the branes and anti-branes wrapping n ho- 
mologous S'^'s of the original geometry are replaced by fluxes threading the n topo- 
logically distinct S'^'s of the new geometry. The local Calabi-Yau threefold after 
the transition is defined by equation ([2|) where the coefficients of correspond 
to the n normalizable complex deformation parameters of the Calabi-Yau. This 
complex equation defines a two-sheeted Riemann surface fibered over the u and v 
coordinates. See figure [3] for a depiction of this geometry. 

The n S'^'s correspond to n 3-cycles Ai such that Ai Ci Aj = for all Dual 
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Figure 3: Depiction of the complex x-plane corresponding to the Riemann surface 
defined by equation ([2]) with uv = 0. The compact A-cycles reduce to counterclock- 
wise contours which encircle each of the n branch cuts of the Riemann surface. The 
non-compact -B-cycles reduce to contours which extend from x = Aq on the lower 
sheet (dashed lines) to x = Aq on the upper sheet (solid lines). 

to each A-cycle is a non-compact 5-cycle such that Ai fl Bj = —Bj (1 Ai = 6ij 
and Bi (1 Bj = for all At the level of the Riemann surface, the Ai reduce 

to n distinct counterclockwise contours encircling each of the n branch cuts of the 
Riemann surface and the Bi reduce to contours which extend from the point x = Aq 
on the lower sheet to the point x = Aq on the upper sheet. The IR cutoff defined 
by Aq in the geometry is identified with a UV cutoff in the open string description. 
The periods of the holomorphic three form Q along the cycles Ai and Bi define a 
basis of special coordinates for the complex structure moduli space: 



where JFq denotes the genus zero prepotential. In the absence of fluxes, each Si 
corresponds to the scalar component of a f/(l) A/" = 2 vector multiplet. Once 
branes are introduced, each Si is identified in the open string description with the 
size of a gaugino condensate. Defining the period matrix: 
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to leading order in the expansion, the Kahler metric for the effective field theory 
is ImTij. For future use we also introduce the Yukawa couplings: 



The net three form flux after the system undergoes a geometric transition is: 

H-i = -f^RR + TiiB-f^NS (9) 

where i^RR is the net RR three form field strength, i^NS is the net NS three form 
field strength and tub = Cq + ie'"^ is the type JIB axio-dilaton. Explicitly, 



N, = J H^, a = a, = -J Hs (10) 

Ai Bi 

for all i. Assuming that the fluxes spontaneously break the M = 2 supersymmetry 
of the system, the flux induced effective potential is [6]: 

{ai + Tn'N''")+-^{N^ + --- + N^) (11) 

where we have introduced a constant shift in order to properly compare the tension 
of branes and anti-branes [6]. Near the conifold point the entries of the r matrix 
satisfy: 

S- 

27rzrij = 5ij log , . / ^ - (1 - ^ij) log — |- + fij {Si, Sn) (12) 



iy(2)(a,)A2 -^^^--A 



where 5ij denotes the Kronecker delta, Aj^ = — aj and fij denotes an analytic 
power series in the variables Si. The degree / terms in fij correspond to the / + 1 
loop contribution to the period matrix in the low energy effective theory. 

Letting B denote the set of branes and B the set of anti-branes, the critical points 
of the one loop effective potential satisfy [6]: 

j^B j^B 
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with vacuum energy density: 



»7r 
9ym 



iG-B jeB 

Solving for the ghieball fields yields: 
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2TTia (An) , 
exp ( ) (A^* > 0) 
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(15) 



Hk 
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1 'V, 1 > 


'Ao^ 


\ 1 



/27rm(Ao)\ 

(16) 



where Q denotes an lA^il**^ root of unity. These discrete phase choices label the 
distinct confining vacua of the theory. 

As explained in greater detail in [6], the brane/anti-brane system decays to the 
supersymmetric ground state by nucleating a bubble of lower vacuum energy den- 
sity. The corresponding domain wall solution is given by a D5-brane wrapping a 
supersymmetric 3-chain in the Calabi-Yau threefold. From the perspective of the 
closed string dual, the D5-brane wrapping the corresponding supersymmetric 3-cycle 
separates vacua with different amounts of flux quanta. 

The type IIA T-dual description of the brane/anti-brane system with two min- 
imal size S'^'s was flrst commented on in [18] and studied in detail in [11]. This 
corresponds to a conflguration with D4-branes and anti-D4-branes suspended be- 
tween a pair of NS5-branes deflned by the equations y = ±W'{x). The analogue of 
the IIB closed string description maps to the regime where the D4-branes and anti- 
D4-branes have dissolved as flux in the NS5-branes. This conflguration was also 
lifted to an M-theory description given by an M5-brane wrapping a harmonic rather 
than holomorphic curve. As noted in [11], although the absence of supersymmetry 
would seemingly preclude any direct connection between the effective potentials of 
the IIB and IIA systems due to the difference in the size of the T-dualized circle, 
the observed match in the limit Qs ^ between the two effective potentials provides 
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further evidence that at leading order in the expansion supersymmetry is indeed 
spontaneously broken by the presence of fluxes. 

Due to the fact that a convenient parametrization of the associated M-theory 
curve is not known for n-cut geometries, we will unfortunately not be able to exploit 
the powerful techniques developed in [11] to describe the Qs ^ limit of the all loop 
phase structure of the associated system for even a restricted set of fluxes. 

Before proceeding to a short review of the two cut phase structure, we first 
comment on the possibility of additional corrections to the Kahler potential due to 
a' corrections. This is a subtle point because in order to keep the string modes which 
mediate supersymmetry breaking from decoupling we must keep a' finite. Although 
in the strict infinite volume limit such effects will apparently not alter the form 
of the Kahler potential, for finite Kahler modulus such corrections will introduce 
important corrections [38]. Because we work in the context of non-compact Calabi- 
Yau threefolds, we shall assume for the purposes of this paper that such corrections 
are negligible. 

2.1 Two Cut Phase Structure 

We now review the phase structure of the two cut system [8]. At leading order in 
an expansion of the periods, the two cut period matrix is [39]: 



As shown in [8], there is an intricate phase structure beginning at two loop order. 
Just as a supersymmetric configuration of A^, branes contains Tr(— 1)^ = N1N2 
energetically degenerate confining vacua, at one loop order the confining vacua of 
the brane/anti-brane system also remain energetically degenerate. These vacua 
correspond to the different roots of unity appearing in the one loop critical points 
of equations f|T5|) and f|T6|) . Beginning at two loop order, this degeneracy in the 
energy densities is lifted. In the closed string dual, the confining vacuum of lowest 
energy corresponds to a configuration where the branch cuts align along a common 
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axis in the complex x-plane. When the size of the confinement scale is fixed to a 
finite value, the action for quantum tunneling to the energetically preferred confining 
vacuum scales as a positive power of N so that the corresponding tunneling rate is 
exponentially suppressed. The two loop contribution to the effective potential also 
destabilizes the vacuum once the flux becomes comparable to the value: 



log 



Ao 



(19) 



Once metastability is lost, the branch cuts begin to expand until they nearly collide. 
When this occurs the interpolating 3-cycle Bi — B2 collapses to nearly zero size. The 
contribution due to additional light states can potentially drive the system to two 
qualitatively different endpoints. Whereas D5-branes wrapping Bi — B2 will tend 
to lower the flux of the corresponding conflguration so that the system relaxes back 
to another metastable conflguration, the contribution due to a D3-brane wrapping 
the same collapsing cycle can trigger a transition to a non-Kahler geometry [8]. 



3 Matrix Model Computation 

In order to study the phase structure of more general geometries, we flrst determine 
the two loop contribution to the period matrix in the n-cut geometry. Although 
it is in principle possible to determine the form of this correction by directly evalu- 
ating the period integrals of the closed string dual, we shall instead use the matrix 
model technology developed in [33-35] to reduce the computation to a perturbative 
Feynman diagram analysis. This will have the added beneflt that we will be able to 
isolate individual string exchange processes which contribute to the vacuum energy 
density. 

The genus zero prepotential of the closed string dual geometry is computed to 
all orders by the planar limit of a large auxiliary matrix model with partition 
function: 

= vwm I '* {-^^''''' '*') 

where $ is a holomorphic N x N matrix and the above matrix integral should be 
understood as a contour integral. The prepotential of the n-cut geometry near the 
semi-classical expansion point is given by expanding the eigenvalues of $ about the 
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n critical points of the polynomial W. The usual eigenvalue repulsion term of the 
matrix model causes these eigenvalues to fill the n cuts of the geometry after the 
geometric transition. With Mj eigenvalues sitting at the z*'^ cut of the geometry, 
this matrix model may be recast as an n-matrix model of the form: 

t. 



■ ■ ■ d^nn exp V {<^ii) Wint ($1 



^MM — — / U'i'ii ■ ■ ■ UM/n^exp I / ^——VVi\^ii) —vyintK^ll,---,^nn 

J]Vol(t/(M,)) 

i=l 

(21) 

where $jj denotes the Mj x Mj block of $ along its diagonal. The periods of the 
yl-cycles are given by the partial 't Hooft couplings of the matrix model: 

Si = gsMi. (22) 

Evaluating Zmm in the saddle point approximation, the genus zero prepotential is: 

3 non— pert ~l~ port (^3) 

where J?-non-pert corresponds to contributions from the Vol([/ (Mj)) factors and ^pert 
corresponds to the perturbative contributions from planar Feynman diagrams: 

" 1 S- 

27ri J'non-pert = X] O*^^ "T^ (24) 



i=l 







27rzJ-pert = 5^ - 5iiy(ai) + C,,...,^S\' ■ ■ ■ S'- (25) 

i=l 0<jl,...,jn 

where the Ci^,,,i^ are coefficients which are in principle calculable. 

We now compute all contributions to J^pert proportional to MjMjM^ for all z, j, k. 
Following [36], the gauge- fixed matrix model action is given by the sum of two 
contributions: 

5'mm = + 5'ghost (26) 

where: 

= -72^1^ — H — ^ ^ ^ 

r>0 i 
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Figure 4: Depiction of all Feynman diagrams which contribute to the MfMj term 
of J-'pert for i 7^ j- In each diagram, ghost propagators are denoted by white and $ 
propagators by dashed lines. The index loops of each diagram are also shown. 



and: ^ 

5'ghost = —^Tr {AjiBijCji + Bij^jjCji + Cij^jjBji) (28) 

where Cij denotes a scalar ghost and Bij its conjugate. In the above we have also 
introduced a constant shift in the definition of the $m's. Suppressing all matrix 
indices, the propagators of the fields are: 

{Bi,C,i) = 1^ (30) 

for all i 7^ j. The /-point interaction vertex of / fields has weight — ly(^) (aj) /gs{l — 
1)! and the three point interaction vertex between the B, C and $ fields has weight 
-l/9s. 

For i ^ j, the three Feynman diagrams which contribute to the MfMj term of 
the free energy are depicted in figure HI Proceeding from the right-most diagram to 
the upper left diagram in a clockwise direction, the contribution of each diagram to 
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the matrix model free energy is: 



For i ^ j ^ k ^ the diagrams which contribute to the MiMjMk term of the 
free energy possess the same ghost and $ field content as the upper left diagram of 
figure m but with all three index loops distinct. The net contribution to the matrix 
model free energy from this class of diagrams is: 

= iy(^)(a.)A.,A.. + - ^) + - (32) 

Two of the three diagrams which contribute to the Mf term of the free energy 
are given by replacing the ghost content of the diagrams of figure H] by $ fields and 
interaction vertices. When three or more cuts are present, a planar figure eight 
diagram also contributes. Starting from the right-most diagram of figure H] and 
ending with the figure eight diagram, the net contribution to the matrix model free 
energy from this class of diagrams is: 

24H^(2)(ai)3 8iy(2)(a,)3 l2W^^){aiY ' 

The two loop corrections to the period matrix are summarized in appendix A. 



4 Effective Potentials and Unoccupied /S^'s 

In this brief section we explain how to treat the effective potential for brane/anti- 
brane systems with unoccupied minimal size S^'s. To this end, we first review the 
treatment of similar supersymmetric brane configurations developed in [40] . Due to 
the fact that all of the initial S'^'s in the open string description are homologous, even 
an unwrapped S"^ will naively undergo a geometric transition to an S^. Note, how- 
ever, that because there is no flux to support the counterparts of the unwrapped 
S'^'s, additional light degrees of freedom corresponding to D3-branes wrapping the 
collapsing S'^'s will enter the low energy dynamics. The effective superpotential now 
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includes contributions of the form: 

aS + V2QLQRS dWeff (34) 

where Qi and Qr denote two M = 1 chiral muhiplets which combine to form the 
J\f = 2 hypermultiplet describing the D3-brane. When the Q fields condense, the 
original disappears and is replaced by an S*^ of size (QlQr) = —Oi/^/2. The 
glueball superpotential for such a system is computed by a matrix model with action 
given by W{x) where the S'j's corresponding to zero flux are set to zero in the 
prepotential. 

In the non-supersymmetric brane/anti-brane configurations studied in this paper 
the effective potential is controlled by the period matrix Tij. Due to the fact that J\f = 
2 supersymmetry is spontaneously broken at leading order in 1/N by the presence 
of positive and negative flux, the prescription for supersymmetric configurations 
described above will continue to hold in the non-supersymmetric context. Labelling 
the fluxes so that Ni,...,Ni ^ for some positive integer / < n with all other A^j 
equal to zero, the appropriate period matrix to use in equation ffTTl) is given by the 
/ X I truncation of Tij to 1 < i, j < I where all glueball fields corresponding to zero 
flux have been formally set to zero in the prepotential. 

5 Two Loop Corrections to the Energy Density 

At one loop order, equation (031) implies that the discrete phase choices for the glue- 
ball fields present in equations (ITSll and (fT6l) have identical vacuum energy densities. 
As in the two cut case, we find that two loop effects lift this degeneracy. The de- 
cay to the confining vacuum of lowest energy proceeds via D5-branes wrapping the 
74-cycles in the closed string dual geometry. We refer the reader to [8] for further 
discussion on the alignment of glueball phases due to tunneling processes. In this 
section we study the orientation of the branch cuts in the lowest energy confining 
vacuum as a function of the roots of W'{x). 

In the remainder of this paper we let St denote the two loop correction to the 
period matrix. Expanding Vcs to linear order in 6t and applying equation (fT3|l 
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yields the first order correction to the vacuum energy density due to two loop effects: 



E = + ^ 4: \Nj\ \Nk\ Re (iSTjk) + 0{St^) (35) 

jeB,keB 

where Srjk is evaluated at the one loop corrected value of the Si's given by equations 
(fT5l) and ( |T6l) . In the following we shall refer to the second term as 6E. Because the 
components of 6Tjk are linear in the Si's, we conclude that the two loop contribution 
lifts the degeneracy in energy present at one loop order. Note in particular that the 
energy only receives contributions from off-diagonal brane/anti-brane components of 

5Tjk. 

From the perspective of the closed string dual geometry, the phases of the Si's 
correspond to the direction of alignment for the branch cuts in the complex x-plane. 
Expanding the defining period integrals for the Si in terms of (a^ — a^) yields: 



S. = ^ 
2m 



a. 



n ^ 

JJ(x - a+)(x - ar)dx = -W^^\ai) {5xif + 0{5x^) (36) 



1=1 



where 4 {6xi)'^ = {a^ — ajY. 

Although there is an ambiguity in the sign of the above period integral due to 
the presence of the square root, the different choices of signs correspond to distinct 
orientations of the A- and B- cycles in the Calabi-Yau threefold. Indeed, once the 
orientation of the -B-cycles is fixed, the choice of signs in the above period integrals are 
then completely fixed. To determine the appropriate sign, we consider a geometry 
given by branch cuts which are all located on the real axis of the complex x-plane 
with g, Ao > 0. Ordering the points so that > > ■ ■ ■ > > a~, we have: 



Ao 

9 



Hi 

2m 



^ R 

W{x - at){x - ai)dx = — (37) 
\ i=i 



where i? is a real number. Adhering to the sign conventions used in [39], general 
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monodromy arguments yield: 



We therefore conclude that the proper choice of sign for 5*1 and therefore all of the 
Si's is the one given by equation fl36l) . 

In the rest of this section we study the alignment of branch cuts in the lowest 
energy confining vacuum. Because there are a discrete number of confining vacua, 
a generic change in the parameters will typically induce a small jump of 0{1/N) 
in the preferred confining vacuum. At large such small jumps do not produce a 
dramatic change in the physics because we already approximate the discrete phase 
choice by a continuous variable. Even in this hmit, however, we find that an ap- 
propriate variation of parameters can induce large discrete jumps in the preferred 
confining vacuum. In subsection 15.11 we treat geometries with three minimal size 
S'^'s and present two examples where the alignment of cuts naturally generahzes the 
results found for two cut geometries. We next demonstrate that for more general 
configurations, changing the relative positions of the a^'s can induce discrete jumps in 
the preferred ahgnment direction of the branch cuts. In subsection 15.21 we generahze 
this analysis to configurations with additional minimal size 5'^'s. 

5.1 Alignment with Three Minimal S^''s 

We now study the energetics of branch cut alignment in geometries with: 

W'{x) = g{x - ai){x - a2){x - as). (39) 

Our expectation is that configurations with a high degree of symmetry will exhibit 
behavior which is similar to that of the two cut system. 

As a first example, we take the a^'s to form an equilateral triangle with A^i > 
D5-branes at x = ai = a, A^2 > D5-branes at a; = 02 = —a, and A^3 < anti-D5- 
branes at x = = i\^a. Setting 6xi = rjC*^' with > 0, the two loop contribution 
to the energy is given by the second term of equation fl35|) : 



SE = Re 



™^ (-1 - 7zVs) ry^^^ + ™^ (-1 + 7zVs) ry^^^ 
™ (11 - 3^V3) + ™ (11 + 3^V3)) rle^^^ 



(40) 
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Minimizing 6E with respect to the /5j's, it now follows that in the lowest energy 
configuration the cuts always point towards the interior of the triangle. Bisecting 
each 60° angle of the equilateral triangle, we further find that the cuts threaded by 
positive fiux align on the side of the bisection closest to the cut with negative fiux. 
Finally, it follows from the last line of equation fHOj) that the orientation of the cut 
with negative fiux depends on the relative magnitudes of |A'^i| and |A''2|. 

In fact, the position of minimal size S'^'s not wrapped by branes will also infiuence 
the alignment of the cuts in the lowest energy confining vacuum. To this end, 
consider a configuration with A'^i > D5-branes located at x = Oi = a > 0, < 
anti-D5-branes located at x = 02 = —a and with a vanishingly small number of 
branes or anti-branes wrapping the 5*^ at a; = 03 = c. In this case {6x3)'^ = and 
the correction 6E takes the form: 

Note that in the limit c — >■ 00, 

lirnSE = -i^^^^Re ((Sxif + (5x2^) (42) 
c-»oo Ana'' 

the corresponding S*^ at x = c decouples from the dynamics of the theory and we 
find that just as in the two cut geometry, the branch cuts align with the real axis of 
the complex x-plane in the minimal energy configuration. 

The symmetries of such configurations still constrain the direction of alignment. 
Indeed, when the Oj's form an isosceles triangle with c = iL a pure imaginary number, 
we find: 

6E = MM Re (Qrle''^^ + Qrle^'^A (43) 
vr 

where 5xi = r^e^^^ and Q is a complex number. It thus follows that the system has 
lowest energy when /3i = —(32 so that the branch cuts point symmetrically towards 
X = c. Indeed, expanding equation ( 14T|) for small L/a yields: 

5E = -MlIM Re (e8*V5a^2g2^A ^ ^-SiL/,a^2^2^p,^^ ^ fli_\ (44) 

so that both cuts tip towards x = c with angle AL/ba. Although the cuts tip so as 
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Figure 5: Depiction of branch cut orientation in the lowest energy confining vacuum 
of the three cut system with iVi > units of fiux through the cut near x = a, N2 < 
units of fiux through the cut near x = —a and N3 = units of fiux through the cut 
near x = iL with L, a > 0. The presence of the additional minimal size S"^ at x = iL 
causes the two cuts supported by fiux to tip up. 



to touch one another, they do not align at the proper angle to touch the point x = c. 
This indicates that when all three cuts are of finite size, the energetically preferred 
configuration corresponds to the case where the endpoints of the cuts are maximally 
close to touching. 

The above examples suggest that although the alignment of the branch cuts in 
the minimal energy confining vacuum depends on the roots of W'{x)^ up to small 
discrete jumps of order 1/7V, this alignment is a smooth function of the relative 
position of the Oj's. Indeed, the coefficients multiplying {dxif' are rational functions 
of the Cj. Note, however, that varying the Oj will produce a discrete jump in the 
alignment of the branch cuts whenever crossing through the zero set of this rational 
function. When the position of a single root of is varied, this zero set reduces to a 
finite number of points in the complex x-plane. In this case a jump can be avoided 
by traversing a contour which does not pass through such a root. 

To establish the existence of such jumping phenomena, we restrict to configura- 
tions where all three a^'s lie on the real axis of the complex x-plane and A'^i > > iV2, 
and A^3 is vanishingly small. With notation as in the previous example, the branch 
cuts in the minimal energy confining vacuum align parallel to the real axis when c 
lies within the union of the following intervals: 

c G (-00, -a) U Q (^4 - Vil) , a) U (^^ (^4 + ^41 j ,00) ^ 5xi G M (45) 

cG (^-cx),-^ (^4 + v/41)) U (^-a,^ (^-4 + v/41)) U (a, cx)) ^ 5x2 G M (46) 
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and perpendicular to the real axis for all other values of c G M which do not cause 
the coefficients of {dxif' to vanish. Hence, there exist regions on the real line where 
both cuts are parallel, both are perpendicular, and one is parallel while the other 
is perpendicular to the real axis of the complex z-plane. Due to the fact that for 
5 > 0, W^'^\a) is positive for c < a and negative for c > a we conclude that the 
glueball phases are also sensitive to the location of the unoccupied S'^. 

Although the locations in the complex x-plane where cut alignment phase tran- 
sitions occur appear to have little geometric content, the individual string exchange 
diagrams determined by the matrix model reflect the underlying location of the D5- 
branes and anti-D5-branes in the open string description. Without loss of generality, 
we focus on the contributions to the alignment of the cut near x = ai = a. Pro- 
ceeding from the right-most diagram to the upper lefthand diagram in a clockwise 
direction in flgure HI the sign of the contribution of each diagram to the alignment 
of 6xi along the real axis of the complex a;-plane is: 



sign 
sign 
sign 



iy(^)(ai) 1 , 

1^(2) (ai)2 Ai2 



(-^7^^2l^W('\a,))^6x^eR forc^K,ai + Ai2] (4 
V W^^> (ai) ^12 / 



( ^^'^^""'O ^^^^^^ ^""^ ^ ""^i ■ ^^^^ 



Note in particular that for the second 1 — 2 exchange diagram the alignment of the 
cut depends on the location of an "image charge" at x = ai + A12. There is a flnal 
class of 2 — 3 string exchange diagrams which can also contribute to the alignment 
of 6x1 when > 0. These correspond to diagrams with the same shape as in the 
upper lefthand diagram of flgure HJ but with all index loops distinct. In this case 
all three diagrams contribute with the same sign to the alignment of 6x1 along the 
real axis: 

sign c > ai. (50) 

VVK^^^ (ai) A12A13 J 

The above analysis indicates that the sign of each string exchange diagram depends 
only on whether c falls to the left or right of x = 02 and x = ai. Indeed, the classical 
trajectories of the string exchange processes in the complex x-plane sometimes must 
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bend in order to avoid passing too close to the S"^ at x = c. 



5.2 Alignment with Multiple Minimal ^^'s 

We now extend the analysis of the previous subsection to configurations where W'{x) 
has multiple roots. Although it is tempting to speculate based on the examples of the 
previous subsection that a large discrete jump in the preferred alignment direction 
of a branch cut will not occur in highly symmetric configurations, in this subsection 
we show that even when the roots of W'{x) appear symmetrically along the real axis 
of the complex x-plane, the orientation of the cuts in the lowest energy configuration 
still suffers discrete jumps. 
Most generally, we have: 



W'{x) = gix' - a') fix) 



(51) 



where f{x) is a monic polynomial with isolated roots. With A^^i > D5-branes 
located at x = a > and A^2 < anti-D5-branes at x = —a with a vanishingly small 
number at the remaining roots of W'{x), the two loop contribution to the energy 
density is: 



A-KO? 



( 



\ 



-3 - 2TTa^ - y ^ I (5. 



i=l 

m 



i=l 
m 



Xx, 



_3 _ 2TTa+c, _ a 



X2) 



i=l 



i=l 



\ 



J 



(52) 



where ci, Cm denote the roots of f{x). Treating the q as m complex coordinates in 
C"*, note that the coefficients multiplying {6xi)^ and (6x2)'^ both vanish on complex 
manifolds of dimension m — 1. It follows that crossing this zero set will induce a 
jump in the branch cut orientation of the energetically preferred confining vacuum. 

Although an analysis indicating the various phase regions for cut alignment for 
general Cj will likely be complicated, when a higher degree of symmetry is present, 
the analysis remains tractable. Restricting to f{x) an even or odd polynomial so 
that: 



X 



X 



(53) 
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with bi all real and a = or 1, we find: 




= Re I I -5 - V-— — -I I {{Sx,r + {6x,f) ] . (54) 



The alignment of cuts is parallel or perpendicular to the real axis of the complex 
x-plane in the energetically preferred confining vacuum when: 



.1 f-rT^>^ (parallel) (55) 

1 1 - [bi/a) 



2 < (perpendicular). (56) 

(bi/a) 



Note that a perpendicular alignment of the cuts requires the existence of a root of 
/(x) such that \bi\ > \a\. This suggests that in addition to the Coulomb attraction 
due to brane/anti-brane forces, the geometry itself exerts a "torque" on the cuts. It 
is tempting to speculate that one possible source of torque may be the condensation 
of D3-branes wrapping collapsing S'^'s which are not threaded by flux. 

As a final example, we consider the alignment of two branch cuts where the 
background geometry is determined by the relation: 

TTX 

W'{x)=gsm—. (57) 

This corresponds to a geometry with an infinite one-dimensional lattice of minimal 
size S'^'s located at x = jL for all j G Z. For configurations with A^^i D5-branes at 
X = mL and N2 anti-D5-branes at x = 0, the two loop contribution to the vacuum 
energy density is now: 



SE 



Re ((-1 + 2(-l)™) {i6x,r + i6x,r)) . (5^ 



Thus, for m even (resp. odd) the cuts align perpendicular (resp. parallel) to the real 
axis of the complex x-plane. 
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6 Geometry and Strong CP 



Upon embedding our configuration in a compact Calabi-Yau, the effective ^-angle of 
tlie brane/anti-brane system becomes a dynamical field with potential determined 
by the vacuum energy density of the open string system. In this section we assume 
that the time-scale for the fluctuations of the corresponding axion is sufficiently 
long compared to the time-scale associated with the decay of a glueball phase to 
the energetically preferred confining vacuum. The axion potential is then given by 
minimizing the vacuum energy density over all of the confining vacua of the theory. 
The roots of W'{x) determine the level of strong CP violation in each strongly coupled 
gauge group. At one loop order the effective 6'-angle of the i^^ D5-brane is defined 
by the relation: 

S, = A',=Q\A,\'exp{te,/m) (59) 

where Q denotes an |A^i|"^ root of unity. The effective ^-angles now follow from 
equations (fT5|) and (fT6|) : 

= m e.TgW^'\a,) - l^.l arg A^. + ^ \Nk\ arg A^, + A {N, > 0) (60) 
e, = maTgW^^\a,)- J2 l^j|argA2, + ^|iV,|argAj-A (iV, < 0) (61) 

where in the above we have absorbed all Aq dependence into the single variable A: 

A = Y,Ni arg A^ + ^ym ( Aq) • (62) 

i 

Note that upon rephasing the parameters W'^'^\ai) and A^^, the value of 9i also shifts. 
This is consistent with the presence of the axial anomaly. The discrete symmetry 
CP is preserved by the z*^ strongly coupled gauge group when 9i vanishes. 

Minimizing the vacuum energy density over all of the confining vacua yields the 
axion potential for the effective 6'-angles: 

Kx = + min {5E) (63) 

where the min denotes minimization over all of the discrete glueball phases and as 
before, 5E denotes the second term of equation (153]) . From the perspective of the 
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closed string dual, it is natural to view V^x as a function of n independent variables 
6*1, ...,6n- On the other hand, from the perspective of the open string description, 
there is a single 5*^. It is therefore also natural to treat the parameter A defined in 
equation fl62|) as the only dynamical axion of the theory. In subsections 16.11 and 16.21 
we determine constraints on the minima of by further studying the behavior of 
6E as a function of the 6'i's as well as the Oj's. 

6.1 CP Invariant Submanifolds 

Viewing the locations of the minimal size S'^'s as an n-component vector ~a = 
(ai,...,a„) G C", we now argue that when / of the S^'s are occupied by a collec- 
tion of branes and anti-branes and n — / is sufficiently large, there exists a real 
submanifold of points in C" with real codimension / such that the corresponding 
axion potential achieves a minimum at di = for all i. It follows from equation 
([35D and the exphcit expressions of appendix A that the two loop contribution to 
the vacuum energy density from terms linear in Si is schematically of the form: 



where the Ci{a , N) are rational functions of the Oj. 6E has a minimum at = 
when Ci < 0. The requirement that all branes and anti-branes preserve strong CP 
imposes / real conditions on the parameter space of the a^'s which will generically 
be satisfied provided n — / is sufficiently large. Hence, there is a manifold of points 
of real dimension 2n — I contained in C"^ such that all branes preserve strong CP. 
Note in particular that this result is independent of whether we treat all of the 6'j's 
or simply A as dynamical fields. 

6.2 Discrete Symmetries and Strong CP 

We now show that discrete symmetries of the Calabi-Yau translate into constraints on 
the amount of strong CP violation. Letting a denote the action of the permutation 
symmetry on the labelling of the branes, we consider brane configurations with W'{x) 
a polynomial with real coefficients such that implies Ni = Nj. 

Returning to equation (|35l) and the explicit expression for 6t given in appendix 
A, the coefficients multiplying 5*^ are complex conjugates of the coefficients in the 




(64) 



i occupied 
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analogous contribution from Sa{i)- To establish this, we consider without loss of 
generality the 6'- angle for A'j > D5-branes at x = a^. Note that: 



Q^ = l. i^^i i^^-i [-^ - ^-W^) + ^W^) ) = 

i^s£B,t£B 

The net contribution from terms linear in Si and ^^{j) is therefore: 

6E{e,) = Re{QS,) (67) 
SE{e^^,)) =Re{ClS^(i)) (68) 

with Ci a complex number. 

We now study the minima of the axion potential first in the case where the 6'j's 
are all dynamical fields, and then in the case where only A is a dynamical field. 
It follows from equations (l67Il and (l68l) that in the first case the minimum energy 
configuration satisfies {9i) = —{Oo-{i))- In particular, when is purely real, the 
effective ^-angle vanishes. Next treat \4x as a function of the single field A. In this 
case, the symmetries of the system in tandem with equation fl6Ql) yield the constraint: 

9,-A = A- 9,^) (69) 

so that A = 9i when aj is purely real. The combined contribution to the energy 
density from 9i and 6'cr(i) now follows from equations fl671) and fl68|) : 



5E {9i) + 6E (^.(i)) = \Q\ \S,\ (^cos (^7. + 



A + 27iki\ f A + 2Tik„ii) 

cos -7i + 



(70) 

where in the above, 

7i = argCi + argS*, - (71) 

Minimizing with respect to ki and A;(7(i), the corresponding axion potential has a 
minimum at A = 0. 
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7 Critical Points and Unoccupied /S^'s 



To further study the phase structure of brane/anti-brane configurations wrapping 
S^'s in the presence of additional minimal size S'^'s, we first determine the behavior 
of the critical points of the corresponding effective potential. As shown in [8], 
although it is difficult to directly solve for the critical points of the n-cut effective 
potential Vcs, it is mathematically simpler to solve for the fluxes as a function of the 
critical points of Ves. In [8] these "attractor-like equations" were used to establish 
the gross features of multi-cut geometries. It follows from section H] that the two 
cut attractor-like equations extend to geometries where only two minimal size S^'s 
are wrapped by branes and anti-branes. We now review the two cut attractor-like 
equations and then apply them in a special class of geometries and flux configurations 
where the analysis of the critical points remains tractable. These results are then 
used in section [H] to show that the breakdown in metastability is similar to the two 
cut geometries studied in [8]. 

The critical points of the two cut system satisfy two two-component vector rela- 
tions [81: 
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(72) 
(73) 



where r is the 2x2 period matrix, C is a constant non-zero complex number, and 
we have introduced: 



Pv 
Pw 



± a/c?! — Adid2 

2d2 

d^ ± a/c?! — Adid2 
2di 



(74) 
(75) 



with: 

di = det 



^ 111 ^ 112 
^112 -^122 



d2 = det 



^ 112 ^ 122 
-^122 -^222 



c/t = det 



-^111 -^112 
-^122 -^222 



(76) 
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The ± signs of equations (j74l) and (j75l) are correlated. The requirement that > 
leads to an unambiguous assignment of brane type for each branch. Switching from 
the + to the — branch of equations (17111 and (1751) changes all branes (anti-branes) 
into anti-branes (branes). 

Next consider more general geometries with: 



W'{x) = g{x' - a') fix) 



(77) 



where g,a > and f{x) is an even or odd polynomial whose roots are all isolated and 
real. This corresponds to a geometry where the minimal size S'^'s are situated along 
the real axis of the complex x-plane and are symmetrically arranged with respect to 
the imaginary axis. We take A'^i > branes at x = ai = a and < anti-branes 
at X = 02 = —a with the remaining minimal size S'^'s left unwrapped by branes. 

We now show that when A'^i = — the corresponding effective potential admits 
critical points such that ti = t2- The analysis is very similar to that given in [8] and 
we will therefore only summarize the various steps in the computation. It follows 
from the Z2 symmetry of the configuration that there exists a finite neighborhood 
around ti, t2 = such that for ti = t2 we have: 
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where M = log(AQ/A^). This implies di = d2, 
Returning to equations fl72|) and fl73l) . ti = t2 implies Ni 



dz and thus, 
—N2 and: 



a 



1 + Pn 



Tn + T12 



M - M\ T22 

^ 22_ 



2ni 



(81) 



8 Breakdown of Metastability Revisited 

In two cut geometries with A^^i = —N2 and both branch cuts aligned along the real 
axis of the complex x-plane, there is a critical amount of flux beyond which the 
effective potential ceases to admit critical points for ti small [8]. Viewing the ratio 
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a/Ni as a function of the real modulus t = ti = ^2 > 0, a critical point of this 
function corresponds to the merging of a local minimum and maximum of \4fr. In 
this section we study the analogous situation when additional minimal size S^'s are 
present in the open string description of the system. 

Based on the condition for a breakdown in metastability in two cut geometries 
given by equation (fT9l) . our expectation is that two loop effects will destabilize the 
metastable vacua when the one loop vacuum energy density becomes comparable to 
the energy density of a supersymmetric brane configuration. Testing this expectation 
requires determining the attractor-like equations for more general flux configurations 
and then determining when the effective potential experiences a merger of two critical 
points. Although for general configurations this appears to be intractable with 
present techniques, when the branes and anti-branes of the open string description 
cluster into two groups so that the system is characterized by the brane/anti-brane 
distance Ascp and the internal separation scale within a cluster Ajnter with |Aintcr| ^ 
the analysis reduces to that of the two cut system. Indeed, it follows from 
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general renormalization group arguments that at energy scales above the mass scale 
set by Ainter, each cluster of branes will behave as a single stack of "unHiggsed" 
branes. In this case the expected breakdown in metastability is given by the general 
bound for the two cut system: 



4^(|A^i| + |iV2|) 
9ym 



TV 



\Ni\ |A^2|log 



Ao 



A 



scp 



>^\N, + N,\ 
9ym 



(82) 



where A^i > denotes the net number of branes from the first cluster and A''2 < 
denotes the net number of anti-branes from the second cluster. This is indeed consis- 
tent with expectations based on treating the branes in the geometry as a distribution 
of point charges. For more general configurations the precise bound on the one loop 
contribution is more involved because the channel of brane/anti-brane annihilation 
depends on the relative positions of all branes and anti-branes in the geometry. This 
complicates the estimate of when to expect a breakdown in metastability. 

To avoid such subtleties, in the rest of this section we restrict our analysis to 
configurations with A'^i = > D5-branes at x = ai = a > and A^2 = — A^ 
anti-D5-branes at x = a2 = —a in geometries with: 



W'{x) = g{x' - a') fix) 



^3) 
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and f{x) as in equation ( !53l) . It follows from the results of section [7| that in this case 
there exist critical points which satisfy ti = t2 = t. Using the explicit expressions for 
Tij and Tijk given in appendix B, expanding equation (15T]) to lowest order in r > 
with t = e^'^r yields: 



2iTia 
N 



± (1 + 2{-lf'^fe^''^pr) log {re'^) - log 



Ao 
A 



+ O {r) + 0{r^ log r) (84) 



where the variable (3 is defined by equation (11231) in appendix B. Because the results 
of subsection 15.21 establish that the minimal energy confining vacuum corresponds 
to a configuration where the branch cuts align parallel or perpendicular to the real 
axis of the complex x-plane, it is enough to study the critical points of 2Tiia/N as a 
function of t G M. Further, due to the fact that each value of 9ym will be attained 
by some value of r by changing by a small amount, it is enough to demonstrate 
that Re[27ria/A^] has an extremum as a function of r for = or vr. 
At leading order in r, the critical points of Re \2TTia/N] satisfy: 
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where the variable = for = and = 1 for (p = tx. The existence of a 
critical point therefore requires {—lY^'^^^^(3 > 0. Returning to the definition of /? 
given by equation f ll23p in appendix B, note that this inequality is in accord with 
the condition for parallel {"d = 0) and perpendicular (i? = 1) alignment given by 
equations (!55|) and (!56|) . respectively. 

The critical value of the 't Hooft coupling A = Oym-^ which the system is no 
longer metastable is given by evaluating equation at the value determined by 
equation (15^ . Dropping all subleading logarithms, a crude estimate of this value is: 
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8.1 Endpoints and Further Transitions 

As shown in [8] for the two cut system, once metastability is lost, the sizes of the 
cuts will expand in order to facilitate the annihilation of flux lines. This expansion 
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leads to the collapse of a 3-cycle in the geometry, at which point additional light 
degrees of freedom enter the low energy spectrum. Whereas the presence of D5- 
branes wrapping this collapsed 3-cycle tend to relax the system back to a metastable 
geometry of lower flux, the presence of a D3-brane wrapping this 3-cycle can drive 
the system to a non-Kahler geometry. 

Even though the alignment of cuts for metastable vacua is sometimes perpen- 
dicular and sometimes parallel to the real axis of the complex x-plane, in all cases 
that we have studied numerically the mode of instability drives the cuts to align and 
expand along the real axis of the complex x-plane once metastability is lost. 

As opposed to the two cut geometry, however, for more general flux configurations 
a branch cut supported by flux may first collide with another cut supported by flux 
of the same sign. In this case, the presence of D5-branes wrapping the collapsing 3- 
cycle separate vacua in the four dimensional spacetime which are nearly degenerate 
in energy density. Hence, in comparison to the two cut geometry, a D3-brane 
wrapping this same collapsed cycle will play a more pronounced role in determining 
the endpoint of this transition. Note, however, that generic values of ^ym can also 
prevent such collisions. Indeed, whereas ^ym tends to tip branch cuts supported 
by positive and negative RR fluxes so that they will continue to collide, when both 
branch cuts support flux of the same sign, they tip so as to remain parallel. 

9 Modes of Annihilation 

When W'{x) has two roots, the branes and anti-branes eventually annihilate by 
traversing an interpolating minimal size 3-chain. In this section we demonstrate that 
for more general geometries, new modes of brane/anti-brane annihilation appear or 
disappear as the locations of the minimal size S'^'s change. 

This is mathematically similar to the problem of determining the BPS spectrum 
of A/" = 2 gauge theories geometrically engineered in type IIB string theory on 
geometries of the form: 



where P{x) is a polynomial with isolated roots. The spectrum of solitons given by 
D3-branes wrapping supersymmetric 3-cycles which interpolate between the roots of 
P{x) depends strongly on the relative positions of these roots [37]. As a particular 
example, when the roots of P{x) are all located on the real axis of the complex 
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x-plane and ordered as r„ < ... < ri, the BPS soliton connecting r„, to ri decays into 
a collection of n — 1 BPS solitons joining rj to rj_i for 1 < i < n + 1 [37]. 

We are interested in a special limit of the above analysis where P{x) ~ Wlx^. 
In this case the BPS spectrum will be similar to that determined by two dimensional 
JV = (2, 2) Landau- Ginzburg theory [37]. It follows from arguments similar to those 
presented in [37] that the spectrum of nearly BPS D5-branes wrapping 3-cycles will 
develop instabilities as the roots of W'{x) change. 

In the next two subsections we restrict our analysis for illustrative purposes to 
geometries with three minimal size 5'^'s. We determine the wall of marginal stability 
for a D5-brane domain wall solution wrapping a 3-chain connecting two roots of 
W'{x) to decay to a two domain wall process. We next demonstrate that such 
solutions have minimal tension. This implies that for more general geometries, 
semi-classical brane/anti-brane annihilation proceeds via a multi-stage "hopping" 
process. Finally, in subsection 19.31 we indicate some qualitative features of this 
phenomenon. 

9.1 Spectrum of Nearly BPS Domain Walls with Three Min- 
imal aS^'s 

We now determine the nearly BPS spectrum of D5-branes wrapping interpolating 
3-chains in geometries with three minimal size S'^'s. Given a 3-chain interpolating 
between x = and x = aj, the tension of the corresponding domain wall is well- 
approximated by: 

Th,] = - [ \W'ix)\\dx\ (88) 
9 s J 

where denotes a contour in the complex x-plane. For a single nearly BPS domain 
wall, the tension is: 

T,, = -\W{a,)-W{a,)\ (89) 
9s 

corresponding to a path in the complex x-plane such that the phase of W'{x) remains 
constant for all x G 7ij. 

An important caveat to equation ( l89l) is that it assumes the existence of a single 
nearly BPS domain wall. Letting ai = —02 = a > and 03 = c, note that c = 
implies W{a) = W{—a) so that the purported domain wall solution connecting a and 
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Figure 6: Plot of the walls of marginal stability with three minimal size 5'^'s located 
at fli = 1, 02 = ~1 and = R + il. As the point is varied throughout the R — I 
plane, a single BPS domain wall solution connecting two of the a^'s may decay into 
a two stage process. When as is located in region a), the process as — > ai decays. 
For as in region b), the process a2 ai decays and in region c) the process a2 — > as 
decays. In region d) all processes of the form aj are allowed. 

—a would have zero tension. This unphysical result establishes the absence of such a 
single BPS domain wall in this region of configuration space. We now determine the 
wall of marginal stability for a single BPS domain wall to decay into two constituent 
products. The central charge of each candidate BPS domain wall is: 



The wall of marginal stability corresponds to the locus of non-trivial c values such 
that the phases of the central charges for the candidate decay products align. Setting 
c/a = R + il, the defining equation for the wall of marginal stability is: 



See figure [6] for the region of marginal stability in the R-I plane. Assuming that 
the brane/anti-brane configuration spontaneously breaks A/" = 2 supersymmetry, the 
overall shape of this wall of marginal stability should remain correct up to 
corrections. 



% = W (a,) - W (a,) . 



(90) 




(91) 
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9.2 Tension Minimizing Solutions 



When a single nearly BPS domain wall solution ceases to exist, the total tension 
of the resulting decay products is a global minimum of the tension formula given 
by equation flHHl) . To establish this, let x{t) denote a parametrization of the path 
7i2 which runs from x(0) = ai to x(l) = 02- It follows from the Cauchy-Riemann 
equations that the net tension of the D5-brane wrapping the corresponding 3-chain 
is: 



where we have decomposed W into real and imaginary parts so that W = Wi +iW2- 
Equation fl92l) implies that minimum tension paths in the complex x-plane map to 
piecewise straight line paths in the complex ly-plane. 

To demonstrate that the two product process described in the previous subsection 
minimizes T, we first restrict to the case with 03 = iL and L G M. In the complex 
ly-plane the decay products of the domain wall solution connecting a to —a map 
to a straight line connecting W{a) to W{iL) and a further straight line connecting 
W{iL) to W{—a). The values of W{±a) and W{iL) in the complex H^-plane are: 



so that for g > 0, the locations of W (±a) correspond to reflections across the real 
axis of the complex VT-plane. On the other hand, for —\Ims{R = 0)\ < L < 
(i? = 0)1, Re [W{iL)] > Re [Vr(±a)] so that W{iL) lies to the right of W{±a) 
in the complex l^-plane. While the shortest path between the points W{a) and 
W{—a) in the complex W^-plane corresponds to a straight vertical line, this may not 
correspond to a path in the complex x-plane. Indeed, consider any path in the 
complex x-plane connecting x = a to x = —a. Such a path will necessarily attain 
the value x = iy for some y G M. In the complex IV-plane, the path must therefore 
pass through the value W{iy). Comparing the values of W{iy) and W{iL), we find: 




(92) 








(93) 



(94) 



W (iy) - W{iL) 



9_ 
12 



(L - yf [Qa' + 2y' + (L + yf) > 0. 



(95) 
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Hence, when Re [ly(zL)] > Re [iy(±a)], the piecewise straight hne path in the com- 
plex ly-plane will be shortest when y = L. Performing an arbitrary real deformation 
of 03 away from the imaginary axis, it now follows that for all points which lie within 
the region bounded by the marginal stability curves given by equation (|9T|) . the de- 
cay products of the single nearly BPS domain wall solution minimize the tension 
functional of equation ( IHHi) . 

9.3 Hopping Effects 

Upon wrapping branes and anti-branes on two minimal size 5'^'s, the results of 
the previous subsection show that the additional S'^'s present in the geometry can 
sometimes obstruct direct brane/anti-brane annihilation. In a semi- classical n- 
domain wall decay process, a collection of branes tunnels from one minimal size S'^ 
to the next until finally annihilating with an anti-brane at the n}^ stage. In the closed 
string dual description, flux lines hop from cut to cut before annihilating. The system 
lowers its energy during each intermediate stage due to the one loop contribution 
to the vacuum energy density from terms proportional to log|abranc — c^anti-branc I • 
Because the tensions of the domain wall solutions are all comparable, the rate of decay 
due to hopping is much slower than direct annihilation. See figure [7| for a depiction 
of this phenomenon. In this subsection we discuss some further consequences which 
follow from hopping effects. 

Letting Fj^j+i denote the decay rate per unit volume for branes to tunnel from 
X = tti to X = Cj+i, the rate for an n-stage annihilation event which begins with 
i = 1 and ends with i = n + 1 is given by summing over the inverse rates of each 
individual link in the decay chair0: 

n 

Ket = Erri.+i (96) 

i=l 

SO that a single hopping event can significantly retard the time-scale for a drop in 
the total number of branes in the system. 

Recall that in the dual closed string description there is a critical amount of 

^In this formula we assume that each successive bubble of vacuum can nucleate anywhere within 
the four dimensional Minkowski spacetime. This provides an adequate approximation of the 
decay of the false vacuum because the decay rate of each term in the decay chain is exponentially 
suppressed. 
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Figure 7: The positions of the roots of W'{x) can sometimes obstruct flux hne 
annihilation in the closed string dual geometry. Instead, the flux must perform 
discrete jumps to other cuts before annihilating. In (a) a hopping event where all 
— units of flux transfer to another cut is shown. The effective potential as a 
function of x G M seen by the branes is shown in (b). 
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flux which can thread an before metastabihty is lost. It follows that in certain 
cases a hopping event can trigger a breakdown in metastability when too much flux 
threads an 5*^. Combining this with the discussion in subsection 18.11 this provides 
a completely dynamical mechanism whereby a system can begin in a metastable 
Calabi-Yau geometry and transition to a non-Kahler geometry. 

As a collection of branes hops from one site to the next, the energetically pre- 
ferred direction of alignment for the branch cuts can also change. As a simple 
example, consider the alignment of branch cuts in the geometry given by an infinite 
one-dimensional lattice of evenly spaced minimal size 5'^'s of the type studied in 
subsection 15.21 For a pair of branes and anti-branes separated by m lattice sites, the 
associated branch cuts align parallel to the real axis for m odd, and perpendicular 
to the real axis for m even. As the branes begin to hop towards one another, the 
energetically preferred orientation of the cuts will also change. 

The minima of the axion potential determined by equation fl35|) will also change. 
Indeed, if the number of occupied minimal size S'^'s changes from / to it follows 
from the analysis of subsection 16.11 that when n — I and n — V are both sufficiently 
large, the manifold of points ( ) G which preserve strong CP will change 

from real dimension 2n — I to 2n — I'. 

10 Radial Mode Stabilization and Glueball Phases 

In this section we propose a speculative mechanism which may stabilize the radial 
mode in compact Calabi-Yau threefolds. Interpreting the two loop contribution to 
the vacuum energy density as a radial potential, we now show that when the phases 
of the glueball fields anti-align with the energetically preferred confining vacuum, 
this potential contains a power law repulsion term which can counter the Coulomb 
attraction term present at one loop. 

For simplicity we first consider a two cut configuration with Ni = N D5-branes 
and N2 = —N anti-D5-branes separated by a distance A = ai — 02. Treating ^ym as 
a mode which has already developed an expectation value, we identify the two loop 
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corrected value of the vacuum energy density with the radial potential: 



Kad (A) 



9ym 
20N^ 



An 



A 



log 



An 



A 



-SttVA 



2TTk + 9 



cos ■ 



YM 



27r/ + e- 



N 



cos ■ 



YM 



N 



(97) 
(98) 



where we have introduced the 't Hooft coupling A = Qym-^ ^ ^ integers 
which label the distinct confining vacua of the system. Note that whereas the 
second term produces a logarithmic attraction between the branes and anti-branes, 
the sign of the power law contribution to Vrad (A) crucially depends on the particular 
confining vacuum of the system. Assuming that the time-scale for the fiuctuations 
of the radial mode are sufficiently fast in comparison to the time-seal^ for the decay 
of the glueball phase to the energetically preferred alignment configuration, it follows 
that there exists a stable minimum of V^ad (A) for k, I ~ N/2 given by: 



lA. 



-20 |Ao 



27rk + 6'ym 27r/ + Oyu 

cos — h cos - 



N 



N 



(99) 



Although highly suggestive, it is unclear whether this mechanism is consistent with 
the stabilization of the other moduli of the system. Indeed, it follows from the 
analysis of [8] as well as section [8] that in the two cut geometry the size of the 
glueball fields becomes unstable when equation holds so that: 



log 



Ao 
A 



20U 



1 

20 



A 

a; 



(100) 



where in the last equality we have approximated by its one loop value. This 
expression is of the same order of magnitude as equation fl^IJl) . Assuming k, I ~ A^/2, 



^Whereas in the axion potential discussion we took the scale of fluctuations for 9ym to be 
sufficiently large compared to all such glueball phase alignment decay events, here we assume that 
the radial mode fluctuates more rapidly. As computed in [8] , when the size of the glueball field is 
not exponentially suppressed this decay rate is an exponentially suppressed quantity so that there 
are typically a wide range of energy scales available at which individual modes become relevant for 
the analysis. 
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simultaneously satisfying equation (j99l) and the condition t <t^. requires: 

4 

<2 (101) 

which is not consistent with an expansion in large Aq. Thus, at the level of anal- 
ysis presented here, we cannot definitively conclude one way or the other whether 
misaligned glueball phases will stabilize the radial mode. 

Adding additional ^^'s does not naively resolve the difficulties encountered above. 
For simplicity, we neglect fluctuations corresponding to the motion of unoccupied 
minimal size S'^'s and consider symmetric conflgurations of minimal size S'^'s of the 
type considered in section [HI It follows from equation (185!) that the size of the 
glueball fleld develops an instability at a value comparable to that given by equation 
fllOOp where the factor of 20 is now replaced by |2/3|, with (3 as in equation fll23p of 
appendix B. On the other hand, this same replacement will also occur in equation 

iSD. 

Before concluding this section, we note that strictly speaking, the analysis of 
section O assumes the existence of another mechanism to stabilize the radial mode. 
This is because the axion potential was obtained by minimizing over all glueball 
phases. 

11 Conclusions and Discussion 

Changing the location of minimal size S^'s triggers phase transitions in geometrically 
metastable conflgurations given by branes and anti-branes wrapping homologous 
minimal size S^'s in a non-compact Calabi-Yau threefold. While two loop effects 
appear to generically lift the degeneracy in vacua present in supersymmetric conflning 
vacua, changing the relative positions of the S'^'s causes discrete jumps in which 
conflning vacuum has minimal energy density. The minima of the associated axion 
potential are also constrained by the discrete symmetries of a given conflguration. 
Although the presence of additional S'^'s does not alter the quahtative conditions for a 
breakdown in met ast ability, their presence leads to the existence of walls of marginal 
stability for nearly BPS brane/anti-brane annihilation processes. This generates 
novel phase dynamics where branes "hop" to nearby minimal size S'^'s. Although far 
from conclusive, when the glueball phases anti-align with the energetically preferred 
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vacuum, the two loop correction to the vacuum energy density may stabihze the 
radial mode present in any compact Calabi-Yau threefold geometry. In the rest of 
this section we discuss some implications of this work and some possible avenues of 
future investigation. 

Over suitable time-scales, the ranks of the low energy gauge groups will dynam- 
ically redistribute due to hopping effects. This leads to a statistical distribution of 
gauge groups connected by such events. It is intriguing that for geometries with a 
large number of minimal size S'^'s, hopping would seem to favor gauge groups with 
a large number of low rank gauge group factors. 

Along similar lines, the statistical mechanics of hopping may also be of indepen- 
dent interest. In geometries with a large grid of unoccupied minimal size S^'s, one 
can arrange for a "current" of branes to flow along the grid. It is natural to speculate 
on the existence of suitable conditions leading to superconductivity in this setup. In 
the case of a one dimensional lattice of minimal size S'^'s, the presence of hopping 
also suggests connections with the physics of spin chains, with the unoccupied S'^'s 
playing the role of impurities. It would be interesting to see whether statistical 
correlations in such systems undergo phase transitions as one increases the "doping" 
of the system. 

It may also be important to determine the consequences of hopping from the 
perspective of a four dimensional observer. As one possibility, note that in some 
cases a multi-stage bubble nucleation process will lead to a collision between distinct 
bubbles. While the collision of such bubbles of vacuum has been studied from 
various perspectives, it would be interesting to see how the results of this paper fit 
with these other developments. 

Walls of marginal stability for brane/anti-brane annihilation processes are most 
likely a generic feature of compact Calabi-Yau threefolds, although there are likely 
to be further features specific to the compact case. For example, the corresponding 
domain wall solutions may not factorize as the product of a path on a local Riemann 
surface with an 5*^ fibration. Determining the "mean free path" for such obstructions 
in a compact Calabi-Yau threefold would be interesting. 

There are many additional subtleties present for brane/anti-brane systems on 
compact Calabi-Yau threefolds. In addition to ol corrections to the form of the 
Kahler metric, one must also find a mechanism which stabilizes all of the non- 
normalizable modes of the non-compact case. It is very suggestive that glueball 
fields with phases anti-aligned with the energetically preferred confining vacuum 
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generate a power law repulsion term which can naively stabilize the radial mode of 
a brane/anti-brane system. Because this result is highly contingent on the phase of 
the glueball fields, in the absence of another mechanism, stabilization of the radial 
mode would appear to preferentially select metastable confining vacua with maximal 
rather than minimal energy density. Note that this is somewhat at odds with the 
procedure outlined in section [6] for extracting the axion potential by minimizing over 
all glueball phases. It is likely that further phase structure results from the interplay 
between the stabilization of the axion, the radial mode, and the other closed string 
modes of compact Calabi-Yau threefolds. 



Acknowledgements 

We thank J. Seo for many helpful discussions and collaboration at an early stage 
of this work. In addition, we thank M. Aganagic, C. Beasley, M.R. Douglas, A.L. 
Fitzpatrick, S. Kachru, J. Marsano, K. Papadodimas, M. Shigemori and L.-S. Tseng 
for helpful discussions. CV thanks the CTP at MIT for hospitality during his 
sabbatical leave, where part of this work was completed. The work of the authors 
is supported in part by NSF grants PHY-0244821 and DMS-0244464. The research 
of JJH is also supported by an NSF Graduate Fellowship. 



Appendix A: Two Loop Corrections to Tij 

In this appendix we collect the two loop corrections to derived in section [3l For 
i j we have: 

+ ' ' (-1 - A..^^^ + 2^^^^^ S- (103) 



^ / 4 4 4 \ 

4- (,l^(2)(a,)A,,A,, + iy{2)(a,)A,,A,fc + 1^(2) (a,) A,, A,, ) + 



2\ 



(104) 



39 



and, for terms on the diagonal: 

" log ^(2)(^^)^2 + yy(2)(o^)3 2 1^(2) (a,)2 ' 

As a check on the above computations, we next compare our value of the two 
loop correction to r which we denote by 6t with the known values for the two and 
three cut geometries. When W'{x) = g{x — ai){x — 02), this reduces to: 



which agrees with appendix B of [39]. 

When W'{x) = g{x — ai){x — 02) (x — 03), we finclf]: 

2 / A23 A32\ ^ 

+ ^ ^ -2 + 5-^-2-^ ^2 108 
(7Af2Ai3A23 V A21 A31; ' ^ ^ 

' ^2 + 5^-2^^3 (109) 



^A23Ai2A32 V '"^31 A21 

2m5T!,-^-' = f , (-2 + 5^ - 2^) (110) 

(7A?2Ai3A23 V A21 A31J ^ ^ 

2 / A32 A23\ ^ 

2-^ + 5-^ ^2 (111) 



(7A^iA^3 V A31 A21/ 
A23 A32 



(?Ai2Ai3A23 A21 A3 



^23 \ ^21 ^31 

which agreed with appendix C of [41]. 



^3 (112) 



■^Wc thank ,J. Sco for suggesting this check. 

''in appendix C of [41], the sign of what is referred to as hab should be reversed in the period 
Ila. This is easily checked (and remedied) by appealing to the symmetry of the r matrix entries 
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Appendix B: Two Cut r^j with n Minimal Size /S^'s 

In this appendix we compute the form of the period matrix when: 

W'{x)=g{x-a,){x-a2)f{x) (113) 

where / is a polynomial with isolated roots and all A-cycle periods other than Si 
and 5*2 are formally set to zero, in accord with the discussion in section |H 

Setting ti = SJ{W^^\a,)Al^), it follows from the expressions given in appendix 
A that the relevant matrix entries are: 



2vririi = log ti 



+ 



^^0 



4 + 2IH + 4/^"^^ 



/(ai) /(ai) 



, Ai2/'(ai: 
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(114) 
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(116) 
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When f{x) is an even or odd polynomial whose roots are all isolated and real, the 



under permutation of its indices. 
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entries of Tij and J^ijk simplify further: 



2niTii 27riri2 
27rzri2 27rzr22 



27rijFiii 27r2jFii2 



log ( tl^ ) + 7^1 + /5t2 - log # - i-l)'''^fp (tl + t2) 



A2 



log - (^1 + t2) log ( t2^ j + 7^2 + 

(120) 
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where we have introduced the parameters: 
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